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ABSTRACT 

Context. This is the third of a series of papers devoted to study in detail and with high-angular resolution intermediate-mass molecular 
outflows and their powering sources. 

Aims. The aim of this paper is to study the intermediate-mass YSO IRAS 20050+2720 and its molecular outflow, and put the results 
of this and the previous studied sources in the context of intermediate-mass star formation. 

Methods. We carried out VLA observations of the 7 mm continuum emission, and OVRO observations of the 2.7 mm continuum 
emission, CO (/=l-»0), C'**0 (7=1-^0), and HC3N (i=12-»ll) to map the core towards IRAS 20050+2720. The high-angular res- 
olution of the observations allowed us to derive the properties of the dust emission, the molecular outflow, and the dense protostellar 
envelope. By adding this source to the sample of intermediate-mass protostars with outflows, we compare their properties and evolu- 
tion with those of lower mass counterparts. 

Results. The 2.7 mm continuum emission has been resolved into three sources, labeled OVRO 1, OVRO 2, and OVRO 3. Two of 
them, OVRO 1 and OVRO 2, have also been detected at 7 mm. OVRO 3, which is located close to the C'^0 emission peak, could be 
associated with IRAS 20050+2720. The mass of the sources, estimated from the dust continuum emission, is 6.5 Mg for OVRO 1, 
1.8 Mo for OVRO 2, and 1.3 Mq for OVRO 3. The CO (7=1-^0) emission traces two bipolar outflows within the OVRO field of view, 
a roughly east- west bipolar outflow, labeled A, driven by the intermediate-mass source OVRO 1, and a northeast-southwest bipolar 
outflow, labeled B, probably powered by a YSO engulfed in the circumstellar envelope surrounding OVRO 1. 

Conclusions. The multiplicity of sources observed towards IRAS 20050+2720 appears to be typical of intermediate-mass proto- 
stars, which form in dense clustered environments. In some cases, as for example IRAS 20050+2720, intermediate-mass protostars 
would start forming after a first generation of low-mass stars has completed their main accretion phase. Regarding the proper- 
ties of intermediate-mass protostars and their outflows, they are not significantly different from those of low-mass stars. Although 
intermediate-mass outflows are intrinsically more energetic than those driven by low-mass YSOs, when observed with high-angular 
resolution they do not show intrinsically more complex morphologies. Intermediate-mass protostars do not form a homogeneous 
group. Some objects are likely in an earlier evolutionary stage as suggested by the infrared emission and the outflow properties. 

Key words. ISM: individual objects: IRAS 20050+2720 - ISM: jets and outflows - stars: circumstellar matter - stars: formation 



1. Introduction 

Molecular outflow is an ubiquitous phenomenon during the ear- 
liest stages of the formation of stars of all masses and luminosi- 
ties. During the last decades, many efforts have been devoted 
to the study and description of the physical properties of em- 
bedded low-mass protostars and their molecular outflows (e.g., 
Bachiller |1996t Richer et al. 120001 ). In recent years, high-mass 
molecular outflows have also been studied in detail, and many 
surveys have been caiTied out towards massive star-forming re- 
gions to achieve a more accurate picture of their morphology and 
properties (e.g., Shepherd & ChurchweU 1996al |1996bi Zhang 
et al. 120011 Shepherd |2005l l. In recent years, there has been a 
growing interest for the Young Stellar Objects (YSOs) located at 
the lower and the higher end of the stellar mass spectrum. Less 
attention has been paid to intermediate-mass protostars, with 
masses in the range 2 Mq < < 10 Mq, which can power 
energetic molecular outflows. In fact, only few deeply embed- 
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ded intermediate-mass protostars are known to date, and only a 
very limited number of their outflows have been studied at high- 
angular resolution. 

In order to study in detail intermediate-mass protostars 
and the outflows that they are powering and to compare their 
morphology and evolution with those of low-mass stars, we 
have started a project to carry out interferometric observations 
of dust and gas towards intermediate-mass star forming re- 
gions. The first two regions observed were IRAS 21391+5802 
(Beltra n et al . |2002l hereafter Paper I) and LI 206 (Beltran 
et al. |2006bl hereafter Paper II). In this third paper, we 
present a thorough interferometric study of the intermediate- 
mass protostar IRAS 20050+2720 and of its molecular out- 
flow. IRAS 20050+2720 is a 280 (Froebrich .2005J YSO 
located in the Cygnus rift at a distance of 700 pc (Dame & 
Thaddeus 1985). Distances in the Cygnus region are highly un- 
certain for stellar objects located at Galactic longitudes close to 
90°, as the kinematic distances are not reliable (e.g. Schneider 
et al. |2006] l. However, IRAS 20050+2720 has a Galactic longi- 
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Table 1. Parameters of the observations 
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(a) For the molecular line observations the 1 cr noise is per channel. 



tude of ~66°, being located in a region where the Cygnus rift 
can be distinguished from the Cygnus X complex. Therefore, 
the kinematic distance determined by Dame & Thaddeus ( |1985l l 
from CO line emission is quite reliable. IRAS 20050+2720 is 
surrounded by a cluster of infrared stars (Chen et al. |1997l l. 
Bachiller et al. (1995) mapped a high-velocity multipolar 
CO outflow in the region, previously detected by Wilking et 
al. ( fl989 ). IRAS 20050-H2720 is deeply embedded in a core that 
has been observed in the continuum at millimeter (Wilking et 
al. IT9891 Chen et al. 1997; Choi et al. |1999l Chini et al. UMH 
Furuya et al.[2005 ; Beltran et al. 2006a]l and centimeter (Wilking 
et al. II989I Anglada et al. |I998al l wavelengths, and also in dif- 
ferent high-density tracers, such as HCO^ and HCN (Choi et 
al.|l999), NH3 (MoUnari et al. 1996), CS and N2H+ (Bachiller 
et al. 1995; WilHams & Mvers [T999] l. and '^CO and C'-^O (Ridge 
et al. i2003i The source has been classified as an intermediate- 
mass Clas s O/I p rotostar (Gregersen et al. |1997l Chini et al. l2001l 
Froebrich |2505T l. 

As shown by these previous studies, intermediate-mass out- 
flows appear more collimated and less complex than previously 
thought (e.g. Fuente et al. 1200 II I when observed with high- 
angular resolution. This argues for the need of high-angular ob- 
servations to better understand intermediate-mass protostars and 
their outflows. To accomplish the ultimate goal of our research 
project, which is to study and compare the properties and evo- 
lution of intermediate-mass protostars with those of their lower 
mass counterparts, in this last paper of the series, we put the 
results on the three intermediate-mass protostars studied in the 
context of intermediate-mass star formation. To do this, we have 
compiled information on outflow properties that have been de- 
rived from interferometric observations to better compare them 
with the ones derived for the sources in our study. 

2. Observations 

2.1. OVRO observations 

Millimeter interferometric observations of IRAS 20050-1-2720 
at 2.7 mm were carried out with the Owens Valley Radio 
Observatory (OVRO) Millimeter Array of six 10.4 m telescopes 
in the L (Low) and C (Compact) configurations, on April 30 and 
May 30, 2003, respectively. The data taken in both array config- 
urations were combined, resulting in baselines ranging from 5.6 
to 42.6 kA, which provided sensitivity to spatial structures from 
about 4'.' 8 to 37". The digital correlator was configured to si- 
multaneously observe the continuum emission and some molec- 
ular lines. Details of the observations are given in Table [T] The 
phase center was located at a(J2000) = 20''07'>5'.80, 5(J2000) 
= -i-27°28'58'.'l. Bandpass calibration was achieved by observ- 
ing the quasars 3C84, 3C273, and 3C345. Amplitude and phase 
were calibrated by observations of the nearby quasar J2015+71, 
whose flux density was determined relative to Uranus. The un- 



certainty in the amplitude calibration is estimated to be ~20%. 
The OVRO primary beam is -64" (FWHM) at 115.27 GHz. 
The data were calibrated using the MMA package developed 
for OVRO (Scoville et al. |I993l l. Reduction and analysis of the 
data were carried out using standard procedures in the MIRIAD, 
AIPS, and GILDAS software packages. We subtracted the con- 
tinuum from the line emission directly in the (u, v) domain for 
C' **0 (7= 1 ^0) and HCi N (7= 1 2 ^ 1 1), and in the image domain 
for CO (7=1^0). 

2.2. VLA observations 

The interferometric observations at 7 mm were carried out on 
January 5, 2002 using the Very Large Array (VLA) of the 
National Radio Astronomy Observatory (NRAOU in the D con- 
figuration. The baselines ranged from 5.4 to 140 kA, which pro- 
vided sensitivity to spatial structures from about 1"5 to 37". 
The phase center was set to the position q'(J2000) = 20''07''5'i91, 
5(J2000) = +27°28'58"2. Absolute flux calibration was achieved 
by observing 3C286, with an adopted flux density of 1.45 Jy. 
Amplitude and phase were calibrated by observing 2023+318, 
which has a bootstrapped flux of 0.48 + 0.01 Jy. CLEANed maps 
were made using the task IMAGR of the AIPS software of the 
NRAO, with the ROBUST parameter of Briggs ( 1995) set equal 
to 5, which is close to natural weighting. We used a Gaussian 
taper in the (u, v) domain when making the maps in order to 
enhance the extended emission in the region. The resulting syn- 
thesized beam is 4'.'9 x 3'.'8 at PA =-9°, and the rms noise is 
0.16 mJy beam"'. 



3. Results 

3.1. Dust emission 

The maps of the millimeter continuum emission around the em- 
bedded source IRAS 20050+2720 at 2.7 and 7 mm are shown 
in Fig. [1] The continuum emission at 2.7 mm has been re- 
solved into three different sources, labeled OVRO 1, OVRO 2, 
and OVRO 3. OVRO 1 and OVRO 2 have also been detected 
at 7 mm, but not OVRO 3, for which an upper limit (3 cr) of 
0.48 mJy beam"' has been estimated. OVRO 3 is not so well de- 
fined at 2.7 mm as the other two sources. However, the source 
is located within ~1"6 of the near-infrared (NIR) source 5 of 
Chen et al. ( 119971 see discussion below about the correction to 
the positions given by these authors). The source is probably 
associated with IRAS 20050+2720 (Fig. [TJ, at least at 12 and 
25 jum wavelengths. However, one cannot discard the possibil- 
ity that part of the emission detected at longer wavelengths (60 
and 100 /vm) is associated with OVRO 1 and/or OVRO 2 as sug- 



' NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. Map of the OVRO 2.7 mm (left panel) and the VLA 7 mm (right panel) continuum emission towards IRAS 20050+2720. 
The rms noise of the map, cr, is 1.5 mJy beam"' (left panel) and 0.16 mJy beam"' (right panel). The contour levels are -5, -3, 3, 5, 
7, 9, 11, 15, 20, 25, and 30 times cr (left panel), and -4, -3, 3, 4, 5, 6, 8, 10, and 12 times cr (right panel). The synthesized beams 
are drawn in the bottom left corner. The error ellipse of IRAS 20050+2720 is indicated in black. The dotted circles represent the 
OVRO (left panel) and VLA (right panel) primary beam (50% attenuation level). Black crosses indicate the corrected position (see 
Sect. 13. It of the NIR sources detected by Chen et al. ( 1997 ), and white crosses the position of the 3 mm sources detected by Furuya 
et al. (I2005I I. The white circle marks the position of the 3.6 cm VLA source reported by Anglada et al. ( I1998al) . Black triangles in 
right panel mark the OVRO 2.7 mm continuum position of the sources. 



Table 2. Positions, millimeter flux densities, sizes, and masses of the cores detected towards IRAS 20050+2720 





Position" 


A = 2.7mm 


A = 7mm 


(arc sec) 


Mass'' 
(Mo) 


n{n2f 
(\(f cm"3) 
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(10^^ cm-2) 
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a(nOOO) 
20'" 07™ 


(5(J2000) 
+27° 28' 


(mjy beam"') 


(mJy) 


(mjybeam"') (mJy) 


OVRO 1 


05^90 


59'.' 3 


464 ± 9.3 


65± 13 


2.04 + 0.17 2.88 + 0.17 


3.3 


6.5 


1.3 


1.4 


OVRO 2 


05^90 


47'.' 9 


17.1 ± 3.4 


18.0 + 3.6 


0.97 + 0.17 1.84 + 0.17 


2.3 


1.8 


1.7 


1.1 


OVRO 3 


06'';43 


53'.' 3 


12.9 ± 2.6 


12.9 ± 2.6 


< 0.48^ 


< 6.7« 


1.3'' 


3.0'' 


1.4'' 



(a) Position of the 2.7 mm emission peak. 

(b) Corrected for primary beam response. 

(c) Deconvolved geometrical mean of the major and minor axes of the half-peak intensity contour. 

(d) Circumstellar mass estimated taking into account all the emission inside the 3 cr level, and assuming a dust temperature of 34 K (see Sect. l4.lt . 

(e) Average H2 volume and column density estimated assuming spherical symmetry and a 10% Helium abundance. The densities have been 
computed taking into account the total circumstellar mass and the size of the source at a 3 cr emission level. 

(f) 3 cr upper limit. 

(g) Source not resolved. 

(h) Estimated from the peak intensity. 



gested by the IRAS high-resolution processed (HIRE^ images. 
OVRO 3 is also associated with a C"*0 peak (see Sect.[33]l. All 
this lead us to consider the bump visible at 2.7 mm emission as 
an independent source, that we called OVRO 3. OVRO 1 and 
OVRO 2 are separated -11" (-8000 AU at the distance of the 
region), while OVRO 3 is separated ~9" (-6000 AU) from the 
other two sources. The position, the flux density at 2.7 mm, the 
deconvolved size of the sources (measured as the geometrical 
mean of the major and minor axis of the half -peak intensity con- 
tour around each source), the mass, and the average volume and 
column density are given in Table |2] 

The strongest millimeter source is OVRO 1. This source is 
located at the center of one of the two CO molecular outflows de- 



^ HIRES images are available as single fields (l°-2°) through the 
Infrared Processing and Analysis Center (IPAC). 



tected in the region (see the next section and Fig.|2]i, which sug- 
gests that OVRO 1 is its driving source. OVRO 1 has a core-halo 
structure at millimeter wavelengths, with a core that has a decon- 
volved diameter of ~3'.'3 (-2300 AU) at the half-peak intensity 
contour of the dust emission, and an extended and quite spheri- 
cal halo or envelope that has a size of -1 1'.'3 (-8000 AU) at the 
3 cr contour level. The source OVRO 2 has a deconvolved diam- 
eter of -2'.' 3 (-1600 AU) at the half-peak intensity contour of the 
dust emission. On the other hand, OVRO 3 is not resolved with 
the angular resolution of the 2.7 mm observations. These sizes 
are consistent with the values found for envelopes around low- 
and intermediate-mass protostars (e.g., Hogerheijde et al l 19971 
1999; Looney et al. ,2000, Fuente et al.^01 ; Paper I, Paper II). 
As can be seen in Fig.[Tl there is extended emission suiTounding 
the three sources, which is clearly visible at 2.7 mm but not at 
7 mm. The estimated upper limit (3 cr) for the extended emis- 
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Fig.2. CO (J-1—^0) emission integrated in different velocity intervals, (-7, +1) km s ' for the SHV blueshifted emission (left 
panel; blue contours), (11, 19) km s"' for the redshifted one {left panel; red contours), (-14, -7) km s"' for the IHV blueshifted 
emission {middle panel; blue contours), and (+19, +26) km s"' for the redshifted one {middle panel; red contours), (-19, -14) 
km s"' for the EHV2 blueshifted emission {right panel; blue contours), and (+26, +31) km s"' for the redshifted one {right panel; 
red contours), overlaid on the 2.7 mm continuum emission at a 3 cr level (greyscale). Contour levels are -5, -3, 3, 5, 7, 10, 15, 20, 30 
and 40 times cr, where 1 cr is 0.6 Jy beam"' km s"' for the SHV interval, 0.5 Jy beam"' km s"' for IHV, and 0.25 Jy beam"' km s"' 
for EHV2. The white crosses mark the millimeter continuum positions of OVRO 1, OVRO 2 and OVRO 3. The black dashed lines 
in left panel outline the direction of outflows A, and B (see Sect. 13.2b . The black dotted lines outline the edges of outflow A powered 
by OVRO 1 . The synthesized beam is drawn in the bottom left corner. The dotted circle represents the OVRO primary beam (50% 
attenuation level). 



sion at 7 mm is 0.48 mJy beam '. This extended component has 
a size of -19" (-13000 AU). 

The position of the eight NIR sources detected by Chen et 
al. ( 11997) are indicated as black crosses in Fig. [1] We com- 
pared these positions with those of the 2MASS catalog and 
found that the positions reported by Chen et al. ( 1997) were 
systematically off'set on average by — 0';29 in right ascension, 
and ~1'.'4 in declination. Therefore, before plotting them, we 
corrected the positions. As already mentioned, the only NIR 
source possibly associated with the millimeter sources detected 
with OVRO at 2.7 mm is source 5, whose position nearly co- 
incides with OVRO 3. Chen et al. (1997) also detect a com- 
pact source at 2.7 mm with the Berkeley-Illinois-Maryland 
Association (BIMA) array coincident with OVRO 1 . These au- 
thors estimate a flux density of 30 + 3.5 mJy, which is half 
of that derived with OVRO (see Table |2]i. Although Chen et 
al. ( 119971 locate the millimeter continuum source close to the 
NIR sources 6, 7, and 8, this is not the case after correcting 
the infrared positions. Therefore, this suggests that OVRO 1 and 
OVRO 2 could be deeply embedded sources still undetected at 
NIR wavelengths. The situation could be different for OVRO 3 
because of its possible association with IRAS 20050+2720 and 
NIR source 5. 

The region has also been observed at centimeter wavelengths 
with the VLAby Angladaet al. ( 1998a). As can be seen in Fig.lT] 
the source OVRO 1 coincides within 0'.'5 with a 3.6 cm source, 
with two components separated by ~ 0'.'3 (Anglada et al. 1998a). 

The position of the millimeter source detected by Furuya et 
al. (|20051 l with the OVRO interferometer at 3 mm is indicated 
as a white cross in Fig. [1] The source is clearly coincident with 
OVRO 1 . These authors also find millimeter emission eastward 
of OVRO 1, which could be associated with one of the NIR 
sources. This suggests that another YSO could be embedded in 
the core (Furuya et al. 2005 ). The position of this possible em- 
bedded YSO detected by Furuya et al. ( I2005I I is also marked in 



Fig. [U with a white cross. Our observations do not have enough 
resolution to separate this eastern emission from that of OVRO 1 
and its extended envelope. The total integrated flux density mea- 
sured by Furuya et al. (2005 ) at 3 mm is 25.6 mJy. If we assume 
a dust absorption coefficient proportional to v^, the expected flux 
at 2.7 mm would be ~ 39 mJy, which is about half of the inte- 
grated flux density measured for OVRO 1 at 2.7 mm. This is 
probably due to the fact that the Furuya et al. observations were 
carried out with more extended OVRO configurations (E and H), 
and thus, being somewhat less sensitive to extended emission. 
As a consequence of this, it is possible that part of the extended 
emission surrounding OVRO 1 was filtered out by the interfer- 
ometer, consistent with the negative contours visible in Fig. 1 of 
Furuya et al. ( 20051 

Choi et al. (11999!) observed the region with BIMA at 3.4 mm, 
with a synthesized beam slightly larger than that of our 2.7 mm 
OVRO observations, and resolve the continuum emission into 
four sources, labeled A, B, C, and D. The peak position of source 
A is in agreement with that of OVRO 1. This source is ex- 
tended and elongated towards the south and the southeast, and 
its structure resembles that of the extended emission surround- 
ing OVRO 1, 2, and 3 as mapped by us at 2.7 mm. 

Chini et al. (120011 ) observed the region with the single dish 
James Clerk Maxwell Telescope (JCMT) at 0.45, 0.85, and 
1.3 mm with angular resolutions ranging from ~8 to 15" and de- 
tected three sources in the region, IRAS 20050+2720 MMS 1, 
MMS 2, and MMS 3, together with an extension to the south 
of MMS 1. The source IRAS 20050+2720 MMS 1 is associ- 
ated with OVRO 1, 2 and 3. Beltran et al. (2006a) detect a sim- 
ilar emission with the Swedish-ESO Submillimetre Telescope 
(SEST) at 1.3 mm. The clump number 1 is the one associated 
with OVRO 1, 2, and 3. However, this clump, which has a size 
of 58'.'2, is much larger than the extended 2.7 mm emission or 
MMS 1, as it includes the extended southern emission as well. 
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Fig. 3. Velocity channel maps of the C'**0 (7=1^0) emission (contours) overlaid on the 2.7 mm continuum emission at a 3 cr level 
(greyscale). The systemic velocity of IRAS 20050+2720 is 6 km s"' . The central velocity of each channel is indicated at the upper 
left corner of the panels. The 1 cr noise in one channel is 40 mJy beam"'. Contour levels are -5, -3, 3, 5, 7, 10, 15, 20, and 30 times 
cr. The conversion factor from Jy beam ' to K is 2.10. The synthesized beam is drawn in the bottom left corner of the bottom left 
panel. The white crosses mark the millimeter continuum positions of OVRO 1, OVRO 2, and OVRO 3. The dashed lines outline the 
direction of the CO outflows A and B detected in the region (see Fig.|2]i. 



Table 3. C'^^O (7=1^0), and HC3N (7=12^11) line and physical parameters towards OVRO 1, OVRO 2 and OVRO 3 in 
IRAS 20050+2720 







Vlsr 


AV 


Tb 


J TedV 


6" 
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Position 


Molecule 


(km s-') 


(km s-') 


(K) 


(Kkms-') 


(arcsec) 


(Mo) 


OVRO 1 


C'«0 


4.68 ± 0.04 


1.50 + 0.21 


0.94 + 0.11 


1.50 + 0.07 








HC3N 


4.82 ± 0.02 


1.44 + 0.05 


3.80 + 0.12 


5.84 + 0.18 


9'.'6 


4.3-5.7 


OVRO 2 


C'^O 


6.33 ± 0.04 


1.70 + 0.10 


1.56 + 0.19 


2.83 + 0.12 


7'.'0" 


<4.3-5.8 




HC3N 


5.86 + 0.08 


1.60 + 0.20 


0.85 + 0.13 


1.44 + 0.17 






OVRO 3 


c'«o 


5.61+0.03 


1.92 + 0.06 


2.51+0.17 


5.14 + 0.15 


14 


11-15 




HC3N 


5.41+0.03 


1.30 + 0.10 


2.06 + 0.14 


2.85 + 0.14 







(a) Deconvolved geometrical mean of the major and minor axes of the half-peak intensity contour. 

(b) Estimated from Eq. 5 of Paper II for typical density distributions with p = 2.0-1.5. 

(c) Emission not resolved. Beam size. 



3.2. The molecular outflow: CO emission 

The CO emission towards IRAS 20050+2720 has been pre- 
viously studied through lower angular resolution observations 
by Bachiller et al. (II995I I. These authors discover three bipo- 
lar molecular outflows in the region: an extremely high velocity 
highly collimated outflow oriented nearly east-west (labeled A 
in Fig. 3 of Bachiller et al.); a bipolar outflow of intermediate 
high velocity oriented in the northeast-southwest direction (la- 
beled B); and a standard high velocity bipolar outflow detected 
along the northwest-southeast direction (labeled C). The obser- 



vations reported in the present study significantly improve the 
angular resolution of previous outflow maps, revealing the struc- 
ture of the molecular outflows in detail, as the interferometer fil- 
ters out most of the extended CO emission seen in the Bachiller 
et al. (|1995t low and intermediate velocity maps. 

Figure |2] shows the maps of the integrated CO {J-1—>Q) 
emission in three different blueshifted and redshifted velocity 
intervals. For the sake of comparison, we chose the same veloc- 
ity intervals as Bachiller et al. (1995): a standard-high veloci- 
ties (SHV) interval, with (-7, +1) km s"' for blueshifted emis- 
sion and (11, 19) km s ' for the redshifted one; an intermediate- 
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Fig. 4. Integrated intensity map of the C'^O (7=1^0) (top) and 
the HC3N (7=12— >11) (bottom) emission over the velocity in- 
terval (3,8) km s"' towards IRAS 20050+2720. Contour levels 
are 0.17, 0.34, 0.68, 1.0, 1.5, 2.0, and 2.5 Jy beam-' km s"' for 
both panels. The white crosses show the millimeter continuum 
positions of OVRO 1, OVRO 2, and OVRO 3. The dashed lines 
indicate the direction of the CO outflows A and B. The error 
ellipse of IRAS 20050+2720 is indicated in black. The synthe- 
sized beam is drawn in the bottom right corner of each panel. 



high velocities (IHV) interval, with (-14, -7) km s ' for the 
blueshifted emission and (+19, +26) km s"' for the redshifted 
one; and an extremely-high velocities (EHV2) interval, with 
(-19, -14) kms"' for the blueshifted emission and (+26, 
+31) km s ' for the redshifted one. The other extremely high- 
velocities interval, EHVl, observed by Bachiller et al. d 19951 1 
has not been mapped because our OVRO observations did not 
cover that range of velocities. The systemic velocity, Vlsr, is 
roughly 6 km s ' . As can be seen in this figure, the CO emission 
traces two bipolar outflows in the region, a roughly east-west 
outflow, coincident with outflow A of Bachiller et al. ( |1995l l, and 
a northeast-southwest outflow, coincident with outflow B. The 
outflow labeled C by Bachiller et al. (1 19951 1 and mapped at the 
lowest velocity interval (SHV) has not been detected. As can be 
seen in Fig. 3 of Bachiller et al. d 19951 1, the outflow C, which has 
a PA of about -35°, has very extended lobes, especially for the 
blueshifted emission that reaches distances of >100" from the 
IRAS source. Most of the emission of this outflow lies outside 
of the OVRO primary beam and, therefore, these interferometric 
observations are not sensitive to it. Note that one cannot rule out 
the possibility that part of the central emission of the outflow C 
overlaps that of the outflows A and B, and as it is not possible 
to distinguish any possible emission of the outflow C from that 
of the outflows A and B, we considered the outflow C as not 
detected in this study. 
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Fig.5. Ci**0 (7=1^0) (top), HC3N (7=12^11) (middle), and 
CO (7=1^0) (bottom) spectra obtained at the position of 
OVRO 1, OVRO 2, and OVRO 3 in IRAS 20050+2720. The con- 
tinuum has been subtracted. The conversion factor is 2.10 K/Jy 
beam-' for C'^^O, 2.14 K/Jy beam"' for HC3N, and 3.30 K/Jy 
beam"' for CO. The thick grey profiles in top and middle pan- 
els are the Gaussian fits to the spectra. The dashed vertical line 
indicates the systemic velocity of 6 km s"'. 



3.3. C'**0 emission 

Figure[3]shows the velocity channel maps for the C"*0 (7=1— >0) 
emission around the systemic velocity, Vlsr - 6 km s"', over- 
laid on the continuum emission. The emission integrated over 
the central channels, velocity interval (3, 8) km s"', is shown in 
the top panel of Fig.|4] The emission is elongated in the east-west 
direction and peaks towards OVRO 3. C'^O shows a more ex- 
tended emission than the dust. In particular, it shows an elonga- 
tion in the direction of the molecular outflow A, which suggests 
some interaction between the outflow and the envelope. The 
spectra taken at the position of the continuum peak of OVRO 1, 
OVRO 2, and OVRO 3 are shown in Figure|5] Table[3]lists the fit- 
ted parameters for C'^O (7=1^0) and HC3N (7=12^11). The 
corresponding Gaussian fits are shown in Fig.|5] As can be seen 
in the channel maps and the spectra, the C"^0 emission is associ- 
ated with the three millimeter sources, with the emission towards 
OVRO 1 weaker than towards OVRO 2 and 3. The velocity inte- 
grated flux density of the central emission is ~16 Jy km s ' and 
has a deconvolved size, measured as the geometrical mean of 
the major and minor axes of the half-peak intensity contour of 
the gas, of -14" or -9800 AU. 



3.4. HC3N emission 

Figure |6] shows the velocity channel maps for the 
HC3N (7=12^11) emission towards IRAS 20050+2720, 
around the systemic velocity, Vlsr - 6 km s"', overlaid on the 
continuum emission. The emission integrated over the velocity 
interval (3, 8) km s"' is shown in the bottom panel of Fig. |4] 
HC3N shows a spatial distribution different from that of C'^O 
and the dust. This is probably because the different molecular 
tracers and the dust are tracing material at slightly different 
physical conditions. The C'**0 elongation along outflow A 
is not traced by HC3N, which could indicate that the gas is 
not dense enough to excite the HC3N (7=12-^11) transition. 
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Fig. 6. Velocity channel maps of the HC3N (7=12— emission (contours) overlaid on the 2.7 mm continuum emission at a 3 cr 
level (greyscale). The systemic velocity of IRAS 20050+2720 is 6 km s"' . The central velocity of each channel is indicated at the 
upper left corner of the panels. The 1 cr noise in one channel is 40 mJybeam Contour levels are -5, -3, 3, 5, 7, 10, 15, 20, 30 
and 45 times cr. The conversion factor from Jy beam"' to K is 2.14. The synthesized beam is drawn in the bottom left corner of the 
bottom left panel. The white crosses mark the millimeter continuum positions of OVRO 1, OVRO 2, and OVRO 3. The dashed lines 
outline the direction of the CO outflows A nd B (see Fig.|2]). 



The spectra taken at the position of the continuum peak of 
OVRO 1, OVRO 2, and OVRO 3 are shown in Figure S The 
emission peaks about 4'.'2 east of the continuum peak position 
of OVRO 1, and it is clearly associated with this source. HC3N 
emission is also associated with the other two sources in the 
region, with the emission towards OVRO 2 weaker than towards 
OVRO 1 and 3. Although the emission is quite compact, it 
shows a tail towards the northeast, and it is slightly elongated 
towards the southeast, along the direction of the molecular 
outflow A. The HC3N emission has a velocity integrated flux 
density of ~12.7 Jykms ' and a deconvolved size of 9'.'6 or 
~ 6700 AU. 



fluxes at 450, 850, and 1300 jum. The mass estimated is 6.5 Mq 
for OVRO 1, 1.8 Mo for OVRO 2, and 1.3 Mq for OVRO 3. 
Note that using a lower dust mass opacity coefficient as pro- 
posed by Hildebrand (il983i would result in masses about a fac- 
tor 4 higher, whereas temperatures 10 K lower would increase 
the masses by about a factor 1.5. The circumstellar masses of 
the three deeply embedded sources are similar to those found for 
intermediate-mass YSOs for OVRO 1, and for low-mass proto- 
stars for OVRO 2 and 3. The average H2 volume density, «(H2), 
and column density, A^(H2), of the sources estimated by assum- 
ing spherical symmetry are given in Table |2] 



4. Analysis and discussion of IRAS 20050+2720 

4.1. Mass and density estimates from dust emission 

Assuming that the dust emission is optically thin, we esti- 
mated the masses of the sources (Table |2]) following Eq. 1 of 
Paper II and adopting a dust mass opacity coefficient at 1 1 1 GHz, 
a:iii=0.2 cm^g"' (kq = 1 cm^g"' at 250 GHz: Ossenkopf & 
Henning'1994') and a gas-to-dust ratio of 100. The dust temper- 
ature assumed to derive the masses is 34 K, which is the dust 
temperature obtained by Chini et al. (1200 11 1 by fitting a grey- 
body to the IRAS fluxes at 60 and 100 fim and the SCUBA 



4.2. Tine moiecuiar outfiows 
4.2.1. Morphology 

The outflow A, which has been detected at the three velocity 
intervals, is centered at the position of OVRO 1. It is elon- 
gated in the direction of PA ^115°, and presents a different mor- 
phology for the blueshifted and the redshifted emission. The 
blueshifted emission is more compact than the redshifted one. 
This lack of symmetry between the lobes could be due to the 
presence of OVRO 3. In such a scenario, the redshifted gas 
would flow freely westwards of the powering source OVRO 1, 
while the blueshifted gas, moving eastwards, would interact with 
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the OVRO 3 core that in projection seems to be located in the 
direction of the blue lobe. The CO spectra towards OVRO 3 
in Fig. |5] exhibits a high-velocity blueshifted wing which sup- 
ports the "impact" hypothesis. Such a situation resembles that 
in IRAS 21391 +5802 (Paper I), where the molecular- outflow in- 
teracts with the dense material in the surroundings and gets de- 
flected. The blueshifted emission shows a lobe that shrinks as the 
velocity of the outflow increases, and the direction of the emis- 
sion is roughly the same for all the velocity intervals. On the 
other hand, the redshifted emission is more extended and shows 
a more complex morphology. At standard (SHV) and interme- 
diate (IHV) outflow velocities, the emission seems to trace the 
walls of the cavity excavated by the outflow. As the velocity of 
the outflow increases, the angle of aperture of the conical cavity 
decreases, and the redshifted emission is detected at larger dis- 
tances from OVRO 1 . The farthest redshifted emission is located 
at -35" from OVRO 1 at the IHV interval. CO is detected at the 
center of the lobe only at EHV2. In fact, at EHV2, the redshifted 
emission is not tracing the walls of the cavity. The fact that no 
CO emission is detected inside the cavity for SHV and IHV indi- 
cates that the gas has already been evacuated. This suggests that 
the outflow and therefore the powering source OVRO 1 could be 
in a later Class evolutionary stage. Figure|7]shows the position- 
velocity (PV) plot along the major axis of the outflow A and cen- 
tered at the position of OVRO 1, and the PV plot along the major 
axis of the outflow B and centered at the center of symmetry of 
the outflow (see below). As can be seen in Fig. El for the outflow 
A there is blueshifted emission at negative offsets, and for the 
outflow B redshifted emission at positive offsets. This probably 
indicates that the two outflows are overlapped along the line-of- 
sight near their centers of symmetry where the powering sources 
are embedded or that the outflows are not resolved at the base by 
the current resolution. 

The bipolar outflow B has a PA of ~2Q° and it has been 
clearly detected at SHV. At IHV, only some weak blueshifted 
emission towards the NE is visible. The powering source is 
probably located close to the center of symmetry of the lobes. 
However, although there is some extended 2.7 mm continuum 
emission towards that position, the angular resolution of our ob- 
servations is not high enough to separate the emission of the 
source driving outflow B from that of the circumstellar enve- 
lope surrounding OVRO 1. There are several 2MASS infrared 
sources close to the center of symmetry of the outflow, although 
it is difficult to determine whether one of these could be the 
driving source. The emission of this outflow presents two well- 
defined and elongated lobes, which have an extension of ~30". 

4.2.2. Physical parameters of the CO outflows 

Assuming that the CO emission is in LTE and is optically thin, 
the mass of the gas associated with the outflows detected in 
the region has been calculated by using Eq. 2 of Paper II. The 
high outflow velocities CO is expected to be optically thin (e.g. 
Arce & Goodman ,2001 ). If CO is optically thick in portions of 
the flow, the correction for opacity should be in any case small 
. As the CO spectra are complex and exhibit prominent self- 
absorption in both interferometric (see Fig. |5]l and single-dish 
observations (Bachifler et al. |1995l Zhang et al. 2005 ), it is difli- 
cult to derive the excitation temperature from the CO brightness 
temperature. Therefore, taking into account the assumption of 
LTE, we have assumed that the excitation temperature is equal 
to the kinetic temperature of 29 K derived from NH3 observa- 
tions by Molinari et al. ( 119961 1, a temperature slightly lower than 
the dust temperature derived by fitting the spectral energy dis- 
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Fig. 7. PV plot of the CO (/=1— >0) emission along the major 
axis of the molecular outflow A with PA= 1 15° (left), and of the 
outflow B with PA= 20° (right). The position offset is relative 
to the position of OVRO 1 (left) and relative to the center of 
symmetry of the outflow B, a(J2000) = 20''07*'05!38, 5(J2000) 
= +27°28'59'.'3. Contours are -1.2, -0.9, -0.72, -0.54, -0.36, 
-0.18, 0.18, 0.36, 0.54, 0.72, 0.90, 1.2, 1.8, 2.7, 3.6, 4.5, and 
5.4 Jy beam The horizontal line marks the systemic velocity, 
Vlsr = 6kms"^ 



tribution (SED) by Chini et al. (2001") but similar to the exci- 
tation temperature used for other intermediate-mass molecular 
outflows. We assumed an [H2]/[CO] abundance ratio of 10'* (e.g., 
Scoville et al. 1 19861 1, and a 10% Helium abundance. In Table |4] 
we give the mass, the momentum, kinetic energy, and momen- 
tum rate in the outflows, which have been estimated following 
the derivation of Scoville et al. ( 1986 ). We also report the range 
of outflow velocities for the outflows. Note that the momentum, 
kinetic energy, and momentum rate in the outflow have not been 
corrected for the (unknown) inclination angle, /, of the flow with 
respect to the plane of the sky. In the case of correcting for in- 
clination, the velocities should be divided by sin /, and the linear 
size of the lobes by cos /. It is almost impossible to separate the 
emission of both outflows for the SHV interval, especially close 
to the center of symmetry of the outflows (see previous section). 
Therefore, it is possible that some emission from outflow B has 
been counted as part of outflow A and vice-versa. As already 
mentioned in Paper II, due to some possible filtering of the emis- 
sion by the interferometer, absorption of the emission by the am- 
bient cloud, opacity effects, and the integration range chosen, the 
masses calculated should be considered as lower limits. The val- 
ues of the mass, momentum, kinetic energy, and momentum rate 
are consistent with the values estimated from interferometric ob- 
servations for other intermediate-mass molecular outflows (see 
Table [Sj. The dynamical timescales of the outflows, which have 
been estimated taking the ratio of the maximum separation of the 
outflow lobes and the terminal velocity of the outflow measured 
from the CO line wings, are 3000 and 5000 yr for the blueshifted 
and redshifted lobes of outflow A, and 5000 and 7000 yr for the 
blueshifted and redshifted lobes of outflow B. 
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Table 4. Properties of the CO outflow 





V" 




pc 




^ out 


Lobe 


(km s-') 


(Mo) 


(Mokms-') 


(10^-* erg) 


(Mokm s"' yr"') 


OUTFLOW A 


Blue 


[-19,+1] 


0.10 


0.99 


1.05 


2.8 X 10^* 


Red 


[+11, +31] 


0.10 


1.05 


1.32 


2.5 X 10-* 


Total 




0.20 


2.04 


2.37 


5.3 X 10-* 


OUTFLOW B 


Blue 


[-14,+1] 


0.02 


0.23 


0.24 


6.4 X lO-' 


Red 


[+11, +19] 


0.03 


0.19 


0.14 


4.5 X 10-^ 


Total 




0.05 


0.42 


0.38 


1.1 X 10-* 



(a) Range of outflow velocities. 

(b) Assuming an excitation temperature of 29 K. 

(c) Momenta and kinetic energies are calculated relative to the cloud velocity, which is taken to be Vlsr = 6 km s"'. Velocity not corrected for 
inclination. 



4.3. The dense cores 

4.3.1 . Physical parameters of the dense cores 

As can be seen in the spectra (Fig. |5]l and the fitted parameters 
(Table O, the gas towards OVRO 1 appears to be clearly at a 
different velocity, Vlsr -4.7 kms than OVRO 2 and OVRO 3, 
which have velocities of VLSR-5.5-6kms"'. 

Following Eq. 5 of Paper II, we estimated the virial mass 
of the sources assuming a spherical clump with a power-law 
density distribution p oc r^'', with p - 2.0-1.5, and neglect- 
ing contributions from magnetic field and surface pressure (see 
Table[3]). The virial masses have been computed from C'**0, ex- 
cept for OVRO 1. For this source, it has been computed from 
HC3N as it peaks close to OVRO 1 and is clearly associated 
with it. For OVRO 2, the emission is not resolved, and therefore 
one can only estimate an upper limit. The values estimated for 
OVRO 1, 4.3-5.7 Mq, are slightly lower than the mass derived 
from the continuum, 6.5 Mq. Taking into account that the mass 
of the central protostar has not been considered in this calcu- 
lation, the total mass (circumstellar plus protostar) of OVRO 1 
would be slightly higher than the virial mass. This would be con- 
sistent with the object undergoing collapse, and consistent with 
accretion driving the outflow. For OVRO 2, from the upper lim- 
its computed it is not possible to say whether the virial mass is 
higher or lower than the continuum mass. For OVRO 3, the virial 
mass estimated is clearly higher than the mass derived from the 
continuum. However, it is not possible to compare both masses 
as the C'^O emission, from which the virial mass has been com- 
puted, is much more extended than the continuum emission, 
from which the circumstellar mass has been estimated. 

4.3.2. Envelope mass from gas emission 

As already mentioned in Paper II, the gas abundance relative to 
molecular hydrogen is the main uncertainty of the estimate of 
the mass of the molecular core. Therefore, instead of measur- 
ing the mass of the gas towards the protostars from the C'^O or 
HC3N emission, assuming a given fractional abundance of these 
species, it is better to estimate their fractional abundance from 
a comparison with the mass derived from the dust emission. We 
integrated the gas emission in the same area used to estimate the 
mass of the gas from the continuum dust emission, and assumed 
that the mass towards the continuum sources derived from C'^O 
and HC3N is that derived from the 2.7 mm continuum emission. 

Following Eq. 3 and 4 of Paper II, the estimated 
C'^0 (7=1— >0) abundance relative to molecular hydrogen is 



2.5 X 10~^ towards OVRO 1 and ~ 5 X 10"** towards OVRO 2 
and OVRO 3. Typical C'**0 fractional abundances estimated to- 
wards molecular clouds are 1.7-4x10"^ (Frerking et al. |1982t 
Kulesa et al. '2005V Therefore, even taking into account the un- 
certainties up to a factor of 5 introduced in the mass estimates 
by the different dust opacity laws used, the abundances derived 
towards OVRO 1, 2 and 3 are lower than the typical fractional 
abundances. 

Regarding the HC3N emission, the estimated abundance rel- 
ative to molecular hydrogen is 1 x 10 towards OVRO 1, 
6.5 X 10 ^ towards OVRO 2, and 8.5 x 10 " towards OVRO 3. 
The HC3N abundance found in the literature can vary by up to 3 
orders of magnitude, from ~3 X 10"" (Vanden Bout et al. 1983]) 
to -10"** (de Vicente et al. '2000'; Martin-Pintado et al. 2005l l. 
Therefore, the HC3N abundances derived towards OVRO 1, 
and especially those derived towards OVRO 2 and 3 are in the 
lower end of the range of fractional abundances estimated toward 
molecular clouds. 

5. Intermediate-mass protostars and their outflows 

This present work, which is focused on the properties of 
IRAS 20050+2720 and of its molecular outflow, is the third 
of a series of high-angular resolution studies devoted to 
intermediate-mass protostars and their outflows, including the 
study of IRAS 21391+5802 (Paper I) and IRAS 22272+6358A 
(Paper II). We put the results on IRAS 20050+2720 in the con- 
text of our previous results, and in the context of intermediate- 
mass star formation. We further compare the properties and evo- 
lution of intermediate-mass protostars with those of their lower 
mass counterparts. We have compiled information on outflow 
properties that have been derived from interferometric observa- 
tions in order to be sensitive to the same spatial scales and to 
better compare the properties with those derived for the sources 
in our study. Our sample is therefore not very homogeneous, 
as very few intermediate-mass outflows have been observed at 
high-angular resolution. In any case, we have been able to com- 
pile data of a total of 14 intermediate-mass protostars, including 
the 3 sources studied by us. 

5.1. Multiplicity in intermediate-mass star-forming regions 

The millimeter continuum emission around IRAS 20050+2720 
has been resolved into three millimeter sources (see Fig. [TJ: an 
intermediate-mass source, OVRO 1, and two low-mass sources, 
OVRO 2 and OVRO 3. In addition, the powering source of out- 
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Table 5. Properties of intermediate-mass outflows derived from interferometric observations 
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L 








^dyn 


NIR or 




Outflow 


(kpc) 




(Me) 


(Me) 


(Mo km s"' yr"') 


(10' yr) 


MIR"* 


Ref. 


IRAS 20343+4129 IRSl 


1.4 


1000-3200 


0.8 


0.04 


1.3 X 10"" 


4.2 


Y 


1, 2 


OMCl-S 137-408'' 


0.50 




2 


0.07 


1.2 X 10"^ 


0.11-3 


N 


3, 4, 5 


LI 206 


0.91 


1200 


14.2 


0.09 


7.7 X 10"' 


8.0-20 


Y 


6, 7 


Cepheus E 


0.73 


70 


13.6 


0.13 


2.3 X lO"'* 


4.0-8.0 


Y 


8, 9, 10 


IRAS 21391+5802 


0.75 


235-500 


5.1 


0.14 


1.4 X 10"^ 


2.0-3.5 


Y 


11, 12, 13, 14, 15 


IRAS 20050+2720 A 


0.70 


280 


6.5 


0.20 


5.3 X lO"'* 


6.5-7.3 


7 


16, 17, 18 


OMC-2/3 MM7 


0.45 


76 


0.36-0.72 


0.34 


1.1 X 10"' 


6.3-21 


Y 


19, 20, 21 


IRAS 20293+3952 A 


2.0 


1050'' 


4 


2.0 


2.1 X 10"- 


4.3 


N 


1, 22, 23 


IRAS 22171+5549 


2.4 


1800 


25 


2.7 


1.2 X 10"- 


7 


N 


24, 25 


IRAS 21307+5049 


3.6 


4000 


53 


3.2 


2.9 X 10"' 


13 


N' 


24, 25 


S235 NNW-SSE 


1.8 






4.0 


1.0 X 10"' 


15-20 


N 


26, 27 


S235 NE-SW 


1.8 


1000 


16 


9.0 


1.4 X 10"2 


5.4-14.7 


?■'■ 


26, 27 


IRAS 23385+6053 


4.9 


1500 


61 


10.6 


9.0 X 10"' 


7.3 


N 


28, 29 


HH 288" 


2.0 


500 


6-30 


11 


1.4 X 10"- 


28 


N 


30 



REFERENCES: 1: Sridharan et al. ( 2002J; 2: P alau et al. ( I2007bt ; 3: Genzel & Stutz ki (1989 1; 4: Zapata et al. j2006b : 5: L. A. Zapata (private 
communication); 6: Crampton & Fisher dl974t ; 7: Paper II; 8: Lefloch et al. ( |1996t ; 9: Moro-Martm et al. {200JJ; 10: Noriega-Crespo et al. 
( |2004t ; 1 1: Mat thews ( 119791 1; 1 2: Sara ceno et al. ([1996); 13: de Gregorio- Monsalv o ( |2006l l; 14: Paper I; 15: G etman et al. gODTt ; 1 6: Dam e & 
Thaddeus( fT985Tl; 17: F roebrich ( [2005> ; 18: this work; 19: Takahashi et al. ( |2006] |; 20: Takahas hi et al . ( |2007> ; 21: Nielboc k et al . ( 12003) ; 22: 
Be uther e t al. j2004a| l; 23: P alau et al. ( |2007ab ; 24: Molinari et al. J2002t ; 25 : Font ani et al. ( |2004ab ; 26: Evans & Blair ( fT98T> ; 27: Felli et 
al. ( 12004) ; 28: Mohnari et al. ( fT998) ; 29: Fontani et al. {2004b ); 30: Gueth et al. d^OOT) . 

(a) Y: source detected at A < 8.0 fim by the Infrared Array Camera (IRAC) of the Spitzer Space Telescope, ISOCAM at 7 pm, or by 2MASS; N: 
source not detected at /I < 8.0 pm\ ?: source for which the association between the millimeter and mid-infrared emission is dubious. 

(b) The information was obtained either from the literature of from the Spitzer Science Center Data Archive interface Leopard. 

(c) Because of the difficulty of treating each outflow separately, the parameters correspond to the whole emission, thus including the two outflows 
associated with OMCl-S 137-408 and with HH288. 

(d) The luminosity of IRAS 20293+3952 is 6300 Lq. However, most of its luminosity comes from a Bl star (L~ 5250 L©: Panagia et al. 11973) 
associated with it and not from the intermediate-mass YSO powering the outflow A (Palau et al. 2007^. 

(e) Mid-infrared emission at 7 pm has been detected towards IRAS 21307+5049 by Fontani et al. (I2004bl l. However, the fact that the mid-infrared 
emission peak is offset by ~6" from the millimeter and submillimeter emission peaks led these authors to suggest that the mid-infrared data are 
likely due to the nearby cluster. 

(f) Mid-infrared emission has been observed 1'.'5 to the south of the millimeter peak (Felli et al. 12006) . Such a separation is slightly greater than 
the error of the relative positions, and according to Felli et al. ( I2006t it cannot be firmly established whether both infrared and millimeter sources 
are related. 



flow B has not yet been detected at millimeter wavelengths, or 
its emission has not been resolved from the extended emission 
of the circumstellar envelope surrounding the intermediate-mass 
OVRO 1 (see Sect. 13. 2t . If one also takes into account the mil- 
limeter source detected by Furuya et al. ( I2005I I. located about 5" 
east of OVRO 1 and not resolved by our observations, the emis- 
sion around IRAS 20050+2720 would be associated, at least, 
with 5 YSOs. And, except for OVRO 3 that could be associated 
with IRAS 20050 +2720 and with the NIR source 5 detected by 
Chen et al. ( 119971 see Sect. 13. It . none of them seems to be de- 
tected at infrared wavelengths yet. This situation very closely 
resembles the scenario around the intermediate-mass sources 
IRAS 21391+5802, in the bright-rimmed cloud IC 1396N, and 
IRAS 22272+6358A, in L1206, where the millimeter emission 
has been resolved into an intermediate-mass source surrounded 
by less massive and smaller objects (Paper I; Beltran et al. |2004al 
Fuente et al. 2007; Neri et al. 120071 Paper II). This suggests that 
intermediate-mass stars form in dense clustered environments. 
In other words, the formation of an isolated intermediate-mass 
star if existing appears rare. This is different from the low-mass 
case, where examples of protostars formed in relative isolation 
are known (e.g. LI 157: Beltran et al. 2004 b, HH 211: Palau et 
al.|2006a, or Bok globules: Laundhart et al. |1997l l. The crowded 
environments would explain the complexity of the immediate 
vicinity of intermediate-mass protostars. 

Intermediate-mass protostars are associated with small pro- 
toclusters of low-mass millimeter sources, some of them not 



yet detected at shorter wavelengths. In many cases, they are 
also associated with young and embedded clusters of more 
evolved NIR sources (e.g. IRAS 21391+5802: Nis ini et a L[200T] 
Getm an et al. |2007l L1206: Ressler & Sh ure [19911 Kumar 
et al. [ 20061 IRAS 20050+2 720: C hen et a. [T997] Gute rmuth 
et al. [2005] Kumar et al. [20061 S235: Felli et al. fT997l 
IRAS 22172+5549: Fontani et al. [2004a Kumar et al. 20061 
IRAS 21307+5049: Fontani et al. l2004a. IRAS 20293+3952 
A: Palau et al. 2007a). Interestingly, in some cases these 
NIR clusters contain a large number of more evolved low- 
mass Class II and C lass III sources (e.g. IRAS 21391 +5802: 
Getman et al. [20071 IRAS 20050+2720: Chen et al. [T9971 
IRAS 20293+3952 A: Palau et al. |2007a[ and possibly L1206: 
M. S. N. Kumar, private communication). This means that the 
intermediate-mass protostars and these more evolved low-mass 
sources cannot be coeval. All this suggests that, at least in some 
regions, intermediate-mass protostars would start forming after 
the first generation of low-mass stars has completed their main 
accretion phase. Such a situation appears to be typical of high- 
mass star-forming regions, where embedded massive stars that 
are still accreting material are surrounded by clusters of less em- 
bedded and more evolved low-mass Class I or Class II sources 
(Kumar et aL l2006l l. 

The presence of multiple sources in intermediate-mass star- 
forming regions also explains the complex appearance of the 
molecular outflows in these regions when observed with low- 
angular resolution. This is due to, firstly, the presence usually 
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Fig. 8. Correlations between source and outflow properties for the sample of intermediate-mass protostars and outflows listed in 
Table |5] (a) Outflow momentum flux Font versus bolometric luminosity Lboi- The 'best fit' Fout-^'boi correlation found for Class I 
sources by Bontemps et al. (11996^ is plotted as a dashed line. Dotted lines represent the dispersion of the fit. The solid line presents 
the relation Fom = hyoi/c. (b) Fout versus circumstellar envelope mass Mem- The 'best fit' Fout-A^env correlation found for a low- 
mass Class and Class I sample by Bontemps et al. (|l996l l is plotted as a dashed line. Dotted lines represent the dispersion of the 
fit. (c) Outflow mass Mout versus Lboi- The least square linear fit found for a sample of low- to high-mass sources by Wu et al. C2004] l 
is plotted as a dashed line. Dotted lines represent the dispersion of the fit. (d) Mout versus Menv. The solid line presents the relation 
Mout - Mgnv (e) Font c/Lboi (dimensionless) versus Msnv/L'^oi (Mem and Lboi in units of Mq and Lq). The light and dark grey boxes 
roughly denote the zones where the Class and Class I sources, respectively, are located in Fig. 7 of Bontemps et al. ( 1996). The 
horizontal and oblique lines that join two dots indicate two different values for the same source or outflow property (see Table |5]l. 
Open circles represent sources not detected at /I < 8 yum, filled circles represent sources detected at /I < 8 /im, and stars represent 
sources for which the association between the millimeter and mid-infrared emission is dubious (see Table. [Sj. The three sources 
studied by us in Paper I, Paper II and this work are labeled. 



of more than one outflow in the region (e.g. HH 288: Gueth 
et al. I200n IRAS 21391H-5802: Paper I; IRAS 20050+2720: 
this paper), and secondly, to the stronger interaction of the out- 
flow with the high-density clumps surrounding the intermediate- 
mass protostar (IRAS 21391+5802: Paper I). However, as shown 
by these studies of intermediate-mass star-forming regions, the 
molecular outflows, which are intrinsically more energetic than 
those driven by low-mass stars, are collimated and less complex 
when observed at high-angular resolution. In fact, they appear 
collimated even at low outflow velocities (Paper I, Paper II), 
similar to the low-mass protostellar flow in HH 211 (Gueth & 
GuiIloteau |T999l Palau et al. l2006bl l. 

5.2. Relationships between protostar properties and outflow 
parameters 

Table |5] shows the properties of intermediate-mass outflows and 
their powering sources derived from interferometric observa- 
tions, which have been compiled from the literature. Table|5]also 
gives information on whether the sources have been detected 
or not at NIR or mid-infrared (MIR) wavelengths. As seen in 
the table, only a very limited number of intermediate-mass out- 
flows have been studied at high-angular resolution. In spite of 
the poor statistics of the sample, we compare the properties of 



intermediate-mass YSOs with those of low-mass protostars. To 
do this, we checked whether intermediate-mass YSOs are con- 
sistent with the correlations between source and outflow prop- 
erties found for a sample of Class and Class I protostars by 
Bontemps et al. ( 1996). In particular, we checked the correlation 
between circumstellar envelope mass Menv and the momentum 
rate in the outflow Fout, the correlation between bolometric lu- 
minosity Lboi and Fg^t, and that between the normalized momen- 
tum Fout c/Lboi and Menv/L[J^^ (see Fig.|8]). Following Bontemps 
et al. ( 119961 ), we used the observed momentum rate in the out- 
flow Fobs (see Table|5]), corrected for inclination and opacity. We 
applied a mean correction factor of 2.9 for the inclination, cor- 
responding to a mean inclination angle with respect to the plane 
of the sky of ~ 32.7°, and a mean correction factor of 3.5 for 
the opacity. Therefore, the corrected momentum flux of the out- 
flow Fout, is Font ~ 2.9 X 3.5 X Fobs ~ 10 X Fobs. For the three 
molecular outflows that had already been corrected for inclina- 
tion angle /, we only applied the opacity mean correction factor 
3.5. In addition, we also compared Menv and the mass of the out- 
flow Mout, and Lboi and Mout. The different symbols in Fig. |8] 
indicates whether the sources have been detected at /I < 8 /im or 
not, or whether the detection is dubious (see Table |5]). 

In Fig.lSjj andlHI?, we compare the relationships between Font 
and Lboi or Menv obtained for the intermediate-mass outflows 
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with those obtained by Bontemps et al. d 19961 1. The 'best fit' 
correlation obtained by these authors is plotted as a dashed line, 
and the dispersion of the fit with dotted lines. As can be seen in 
the plot, intermediate-mass YSOs have, in general, higher Fout 
than low-mass objects. One could argue that this is due to the 
fact that the sources in our sample have been observed with in- 
terferometers while those of Bontemps et al. with single-dish. 
However, for those sources in our sample that have been ob- 
served with single-dish as well, the values of Fout derived are 
either comparable or even higher than those derived with inter- 
ferometers. According to Eq. 5 of Bontemps et al. ( 1 19961 1, Fout is 
related to the accretion rate Mace, the entrainment efficiency, and 
the outflow driving engine efficiency. Therefore, the higher Fout 
values suggest either a higher Mace for intermediate-mass stars, a 
higher entrainment efficiency or a higher outflow driving engine 
efficiency. Calvet et al. (|2004| l have found, by means of optical 
and ultraviolet observations, that the average mass accretion rate 
for a sample of intermediate-mass T Tauri stars, which are the 
evolutionary predecessors of Herbig Ae stars, is a factor ~ 5 
higher than that for their low-mass counterparts. Therefore, it is 
very likely that intermediate-mass YSOs accrete material faster 
This seems a reasonable conjecture as these protostars have to 
gain enough mass before the circumstellar disks are dispersed. 
According to Fuente et al. (120011 1. the dispersal of the circum- 
stellar material and the end of the outflow activity should occur 
in < 10^ yrs. 

As seen in Table |5] and Fig. [8] half of the sources have not 
been detected at A < 8.0 yum. This suggests that these sources, 
represented by open circles in Fig. |8] could be more deeply em- 
bedded in the core, and therefore to be in an earlier evolution- 
ary stage than those detected and represented by filled circles. 
The difference in detection rate at NIR or MIR wavelengths is 
also evident in some source and outflow properties. Fig.[8]shows 
that for a similar Lboi or Mgnv, the sources not detected at NIR 
or MIR wavelengths have higher, in some cases up to ~2 or- 
ders of magnitude higher, Fo^t and Mout than those detected. 
Therefore, this suggests that the sources not detected at NIR or 
MIR wavelengths are more efficient at driving their outflows, 
which are more powerful, than those detected. The same hap- 
pens when comparing Mout for objects with similar Lboi (see 
Fig. |8]:). Those not detected at NIR or MIR wavelengths put 
more mass in the outflows than those detected. In addition, the 
non-detected sources agree better with the correlation between 
Mout and Lboi obtained by Wu et al. ( 12004 ') for a sample of low- 
to high-mass YSOs. Although the radiation pressure of the cen- 
tral source would not be enough to drive the outflow if the pho- 
tons emitted by the central source were scattered only once, as 
shown by the fact that all the outflows have Fo^t above the ra- 
diative momentum flux Fi-ad - Lboi/c (see Fig.jSjj), Lboi of the 
central source and Mout are still correlated, showing the depen- 
dence of the outflow on its driving source. As seen in Fig. |8]i, 
in general, the sources not detected at NIR or MIR wavelengths 
put almost the same amount of mass in the outflow as they have 
in the circumstellar envelope. 

To further investigate these differences, we removed any 
luminosity dependence from both Fout and Mdw. Following 
Bontemps et al. (1996), we plotted the outflow efficiency 
Foutc/Lboi versus M^n/L^'^^ (Fig. [8^), along with the zones 
where the low-mass Class and Class 1 sources are roughly lo- 
cated in Fig. 7 of Bontemps et al. ( 119961 1. This plot shows that, 
in general, the sources not detected at NIR or MIR wavelengths 
have a higher outflow efficiency. In the low-mass regime, there 
is a decline of the outflow efficiency with evolutionary stage 
(Bontemps et al. 1 19961 ). Therefore, although we do not claim 



that the sources not detected at /i < 8.0 i-tm are like low-mass 
Class objects, we propose that these sources could be younger 
than those detected, as additionally suggested by the differences 
in properties such as Fom, Mout, and outflow efficiency. In the 
low-mass regime, the youngest protostars have been classified 
as Class objects, based on criteria such as the ratio of submil- 
limeter to bolometric luminosity, Lsubmin/i'boi > 5 x lO"-', where 
Lsubmin is measured longwards of 350 yum (Andre et al. 1993]). 
However, intermediate-mass star-forming regions are located, in 
general, at farther distances than low-mass ones, and usually 
associated with small clusters. Hence, the luminosity at wave- 
lengths longwards of 350 fj.m will likely have contributions from 
other sources in the cluster, and therefore the criterion used to 
classify the youngest low-mass sources might not be applica- 
ble for intermediate-mass ones. In order to confirm this possible 
evolutionary difference in the sources of our sample, one should 
conduct a detailed source by source study of the SED at high- 
angular resolution, however, this goes beyond the scope of this 
work. 

6. Conclusions 

We studied the dust at 2.7 and 7 mm and the gas at 2.7 mm to- 
wards IRAS 20050-1-2720, an intermediate-mass YSO, with the 
OVRO Millimeter AiTay and the VLA. We also put the results 
on this and previously studied intermediate-mass sources in the 
context of intermediate-mass star formation. 

The 2.7 mm continuum emission has been resolved into three 
sources, OVRO 1, OVRO 2, and OVRO 3. Two of them, OVRO 1 
and OVRO 2, have also been detected at 7 mm. OVRO 3, which 
is located close to the C'^0 emission peak, could be associ- 
ated with IRAS 20050-1-2720. The strongest source at millimeter 
wavelengths is OVRO 1 , which has a deconvolved diameter of 
~3'.'3 (~2300 AU) at the half-peak intensity contour of the dust 
emission. The source is surrounded by an extended and a fairly 
spherical halo or envelope, which has a size of ~11'.'3 (~8000 
AU) at the 3 cr contour level. OVRO 2 has a deconvolved diam- 
eter of ~2'.'3 (~1600 AU) at the half-peak intensity contour of 
the dust emission, and OVRO 3 is not resolved. The mass of the 
sources, estimated from the dust continuum emission, is 6.5 M© 
for OVRO 1, 1.8 Mo for OVRO 2, and 1.3 for OVRO 3. 

The multiplicity of sources observed towards 
IRAS 20050-H2720, as well as towards IRAS 21391-1-5802 
and LI 206, appears to be typical of intermediate-mass star- 
forming regions. Intermediate-mass stars form in dense 
clustered environments, and in most cases they are associated 
with small protoclusters of low-mass millimeter sources, some 
of them not yet detected at shorter wavelengths. In many cases, 
as for example IRAS 20050-H2720, IRAS 21391-1-5802 and 
LI 206, intermediate-mass protostars are also associated with 
young and embedded clusters of more evolved NIR sources that 
sometimes contain a large number of more evolved low-mass 
Class II and Class III sources. This suggests that, at least in 
some regions, intermediate-mass protostars would start forming 
after the first generation of low-mass stars has completed their 
main accretion phase, a situation that appears to be typical of 
high-mass star-forming regions (Kumar et al. 2006). 

The CO (7=1^0) emission towards IRAS 20050+2720 
traces two bipolar outflows within the OVRO field of view, 
a roughly east-west bipolar outflow, coincident with the out- 
flow A of Bachiller et al. (1995) and apparently driven by the 
intermediate-mass source OVRO 1, and a northeast-southwest 
bipolar outflow, labeled outflow B and probably powered by 
a YSO engulfed in the circumstellar envelope surrounding 
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OVRO 1 . The outflow A presents a diflferent morphology for the 
blueshifted and the redshifted emission. The blueshifted emis- 
sion is more compact than the redshifted one. This lack of sym- 
metry between the lobes could be due to the presence of the 
OVRO 3 core. High-angular resolution observations show that, 
in general, intermediate-mass outflows, which are intrinsically 
more energetic than those driven by low-mass YSOs, are not in- 
trinsically more complex than low-mass outflows. Intermediate- 
mass outflows appear collimated, even at low velocities, and 
have properties that do not differ significantly from those of 
low-mass stars. Fout values of intermediate-mass protostars are 
higher than those of low-mass ones. Although one cannot dis- 
card a higher entrainment efficiency or outflow driving engine 
efficiency, this suggests that intermediate-mass YSOs very likely 
accrete material faster than low-mass ones. 

Intermediate-mass protostars do not form a homogeneous 
group but have different properties. The sources not been de- 
tected at A < 8.0 /im clearly have higher Fout, Mom, and outflow 
efficiency than the sources detected at A < 8.0 fim for similar 
Menv and Lboi. This suggests that the sources not detected at NIR 
or MIR wavelengths, which are likely more embedded in the 
core and drive more powerful outflows, could be in an earlier 
evolutionary stage. 
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